Solar neutrino measurements in Super Kamiokande-II 
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The results of the second phase of the Super-Kamiokande solar neutrino measurement are pre- 
sented and compared to the first phase. The solar neutrino flux spectrum and time-variation as 
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well as oscillation results are statistically consistent with the first phase and do not show spectral 
distortion. The time-dependent flux measurement of the combined first and second phases coincides 
with the full period of solar cycle 23 and shows no correlation with solar activity. The measured 8 B 
total flux is (2.38±0.05(stat.)±o'.is( s y s -)) x 10 B cm 2 sec 1 and the day-night difference is found to 
be (—6.3 ± 4.2(stat.) ± 3.7(sys.))%. There is no evidence of systematic tendencies between the first 
and second phases. 



I. INTRODUCTION 

The first phase of the Super-Kamiokande experiment, 
SK-I [1-2], yielded high precision measurements of the so- 
lar neutrino flux. In spite of the loss of numerous photo- 
multipler tubes [PMT] sustained in an accident, SK con- 
tinued to collect data with reduced photo-cathode cover- 
age and a higher energy threshold. Data collection and 
analysis methods had to be revised due to the loss of 
detector sensitivity. Super-Kamiokande 's second phase 
[SK-II] ran from December 2002 to October 2005. 

Throughout this paper, the methods and results of SK- 
II are compared with SK-I and, when differing, are de- 
tailed for SK-II. 



II. SK-II PERFORMANCE 
A. Detector Simulation 

It was determined that the November 12th, 2001 acci- 
dent sustained by the SK detector was caused by a prop- 
agating shock wave initiated by an imploding PMT lo- 
cated at the bottom of the inner detector. Therefore, 
blast shields were installed to protect the PMTs against 
such a chain reaction. These shields are 1.0 cm-thick 
transparent acrylic domes allowing light to pass to the 
PMTs' photo-sensitive surface. This presents an ad- 
ditional medium through which Cherenkov light must 
travel. Reflection and refraction of light on the acrylic 
surface is accounted for in the SK-II GEANT 3 Monte 
Carlo detector simulation. The acrylic shields' trans- 
parency at normal incidence is better than 98% above 
400 nm in wavelength. It is about 86% at 300 nm. 

For light propagation in water, both SK-I and SK- 
II adopt a 3-part model of light attenuation consisting 
of Rayleigh scattering, Mie scattering, and absorption. 
We consider two types of absorption: long wavelength 
(A > 350 nm) and short wavelength (A < 350 nm) ab- 
sorption. In the long wavelength absorption region, we 
utilize an independent model derived using direct mea- 
surements from an integrating chamber absorption meter 
[ICAM] applied to pure water [5] . Scattering coefficients 
and absolute short wavelength absorption are tuned to re- 
produce energy distributions in LINAC calibration data 
(see section II. C for a description of the LINAC data). 
At short wavelengths, the SK-I model varies the absorp- 
tion coefficient to describe the changing SK in-tank wa- 
ter transparency as measured by decay electrons from 
cosmic-ray muons. In SK-II, the best description has 



no short wavelength absorption but increased scattering. 
Figure [T] shows the attenuation model in both short and 
long wavelengths for SK-II. 

In determining the expected solar neutrino flux spec- 
trum for a range of oscillation parameters, SK-II follows 
the method of SK-I: the total 8 B and hep flux values of 
the BP2004 Standard Solar Model (SSM) g] are used 
with the neutrino spectrum based on the /3-delayed a 
spectrum of 8 B decay by Ortiz [5] to calculate the flux of 
a particular energy bin. The uncertainties of the neutrino 
spectrum are taken from Bahcall [Ej. 
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FIG. 1: Wavelength dependence of the water parameter com- 
bined SK-II absorption-reflection coefficient. 



B. Event Reconstruction 

1. Vertex 

The determination of event vertex, direction, and en- 
ergy with the reduced light collection capability of SK- 
II has prompted the development of new reconstruction 
methods. For vertex reconstruction, the efficiency of the 
SK-I standard vertex fit significantly drops at energies 
below the SK-I analysis threshold of 5.0 MeV. With 40% 
photocathode coverage, this corresponds to roughly 25 
PMT signal hits. At 19% coverage in SK-II, 25 hits trans- 
lates to 8 MeV. 

The timing residual in an event is defined as the time 
difference between a PMT's hit time fj and the emission 
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time to (fitted to minimize all timing residuals) minus the 
time it would take Cherenkov light to reach that PMT 
given the event's vertex v in the tank: 



^residual — \pi ^o) P 



/C, 



(2.1) 



where hi is the vector location of the hit PMT and c is 
the group velocity of light in water. 

In SK-I, vertex reconstruction is accomplished by se- 
lecting a limited number of hit PMTs from an event (to 
reduce bias from PMT dark noise and scattered light 
hits) and calculating a goodness relation based on the 
timing residuals of the selected hits and a candidate ver- 
tex v. A systematic grid search of candidate vertices is 
performed until the goodness reaches a maximum value. 
After that, the vertex position is fine-tuned to further 
maximize the goodness. The SK-I reconstruction will 
not attempt a vertex fit for less than 10 hits. 

In contrast, the SK-II reconstruction uses all hits from 
an event to form the timing residuals for determination 
of the vertex position. Bias from PMT dark noise is re- 
duced by constructing a likelihood describing the shape 
of the timing residual distribution from LINAC calibra- 
tion data. This likelihood is then maximized from a ver- 
tex search based not on a grid pattern but from a list of 
vertex candidates calculated from PMT hit combinations 
of 4 hits each. The four-hit combinations each define a 
unique vertex given their timing constraints. Any event 
with four hits or more is reconstructed. 

SK-II also makes use of the SK-I goodness-grid search 
method in its online and initial offline analysis for filter- 
ing background events. The final reconstruction, or the 
standard fit based on the residual likelihood method, is 
the final determination of vertex position and can also 
be seen as a correction for any misreconstructed events 
which survived the filtering process. Figure [2] shows the 
various vertex resolutions for the SK-II vertex reconstruc- 
tion. The uncertainty of the measured solar neutrino rate 
due to systematic shifts in vertex position is estimated to 
be 1.1%. 

Both SK-I and SK-II utilize a fast fit online reconstruc- 
tion method for pre-filtering low energy events. Details 
can be found elsewhere [2]. 
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FIG. 2: Vertex resolution (defined as 68.2% of reconstructed 
events which reconstruct inside a sphere of radius a from the 
correct vertex) of 8 B Monte Carlo events as a function of total 
recoil electron energy. 



by detector coverage but by multiple scattering of elec- 
trons in the tank. The difference between data and MC 
angular resolutions, however, is greater in SK-II due to 
the larger discrepancies in energy scales between data 
and MC. Therefore, we assign an angular resolution sys- 
tematic error of 6.0%. Angular resolution is shown in 
Figure [3] 
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FIG. 3: Directional resolution of Monte Carlo events as a 
function of recoil electron total energy. 



2. Direction 



The direction reconstruction is identical to the SK- 
I method: a likelihood function is used to compare 
Cherenkov ring patterns between data and MC distribu- 
tions. Opening angles between the particle direction and 
reconstructed vertex-to-hit PMT position are scanned us- 
ing grid searches at varying levels of precision. The SK-II 
standard fitter is used to determine the vertex. The ab- 
solute angular resolution (defined as the maximum an- 
gular difference between 68% of the reconstructed and 
true event directions as determined by MC) of SK-II dif- 
fers from SK-I by about 10% and is mostly limited not 



3. Energy 

The reconstruction of event energy is identical to 
that of SK-I with modification specific to SK-II (photo- 
cathode coverage, blast shields, etc.). From the number 
of in-time hit PMTs (coincident within 50 ns) from an 
event (N hit ), various corrections are made. The resulting 
effective hit sum has a consistent value throughout the 
detector for a given event (N cff ). From N eff , we deter- 
mine energy. Refer to [5] for specific information on the 
conversion from N ha to N eB . 
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The -/V eJ -tc-energy conversion function must be modi- 
fied for SK-II due to smaller values of N cff corresponding 
to equivalent energies with larger N cff in SK-I. This is 
done by generating MC events at discrete input energies 
between 5 and 80 MeV, calculating their N cff values, and 
then interpolating the energy function. As with SK-I, 
energy refers to total energy of the event (kinetic energy 
plus electron rest mass). 

Since the corrections in N cff depend on the water trans- 
parency, the reconstructed energy also varies slightly 
with changing water quality. See Figure [4] for N eff 
as a function of time for a given water transparency. 
When calculating energy for data events, the water trans- 
parency value as determined by decay electrons from 
cosmic-ray muons is used as an input parameter. How- 
ever, for MC events, the change in water transparency 
is not simulated due to its relative stability compared to 
SK-I and a calculated, constant value of 101 m is used 
for all MC events in the analysis. 
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FIG. 4: Upper figure shows the time variation of the mea- 
sured water transparency (weighed by the Cherenkov spec- 
trum) during SK-II. Lower figure shows the stability of the 
SK-II energy scale as a function of time. The absence of data 
points in late 2003 is from detector dead time due to an elec- 
tronics upgrade. 



An analytical function of the detector's energy resolu- 
tion can be determined with the same MC events used 
for the ./V eir -to-energy conversion function. The energies 
of the MC events are calculated from the method de- 
scribed above and their fitted Gaussian mean energy and 
corresponding 1 sigma values are plotted for each discrete 
energy. A sigma function a(E) is then fitted to use in a 
normal Gaussian probability density function 



P(E,E') 
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2ttct 



exp 



(E 1 - Ef 



(2.2) 



given by 

a(E) = 0.0536 + 0.5200^ + 0.0458£, (2.3) 

in units of MeV. The SK-I resolution is a = 0.2468 + 
0.1492\/i? + 0.0690£. Both resolutions are shown in 
Figure [5j Equation 2.2 with Equation 2.3 can be used 
to apply the SK resolution when calculating theoretical 
spectra for comparison with SK data. 




where E is the electron's true recoil energy and E' is 
the reconstructed energy. The function a(E) for SK-II is 
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FIG. 5: Energy resolution as a function of total recoil electron 
energy of MC events. The red dashed line is SK-I. 



C. LINAC and 16 N Energy Calibration 

As with SK-I, the primary instrument for energy cali- 
bration in SK-II is an electron linear accelerator [LINAC] . 
Detailed discussions on the LINAC calibration methods 
can be found elsewhere [7]. Electrons are injected into 
the SK tank at various positions (see Figure [6]) at ener- 
gies between 5.8 and 13.4 MeV. After reconstructing the 
energies of LINAC events, these data are compared with 
MC to determine the deviation in energy scales. Vari- 
ous MC parameters are then adjusted to minimize the 
differences (see section II. A). 

The minimum uncertainty in the SK-II absolute energy 
scale using 13.4 and 8.8 MeV LINAC data is calculated 
to be 1.4%. This is in contrast to the SK-I estimated 
value of 0.64%. Figure [7] shows the relative difference of 
reconstructed energies of LINAC data and MC as well as 
their differences in energy resolution. 

16 N is also used as a calibration source in conjunc- 
tion with the LINAC calibration [5]. 16 N is produced by 
lowering a deuterium-tritium neutron generator into the 
tank and initiating the fusion reaction 2 H+ 3 H— > 4 He+n. 
A fraction of these 14.2 MeV neutrons collide with 16 in 
the water to produce 16 N which then decays with a half- 
life of 7.13 seconds. In most cases, the Q- value is shared 
between 6.1 MeV gamma rays and a /3-decay electron. 

16 N decays allow directional studies on the energy scale 
not capable with the unidirectional LINAC beam. At a 
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FIG. 6: The LINAC system at SK. The dotted line repre- 
sents the fiducial volume of the detector and the numbers 1-6 
indicate where LINAC data were taken in SK-II. 



total 16 N decay product energy of 10.4 MeV, observed 
energy at several tank positions is compared with MC- 
simulated energy and the difference is shown to agree 
with those obtained from LINAC data and MC. The 16 N 
energy scale difference is averaged to be ±1.2% compared 
with LINAC's ±1.4%. In addition, the isotropic 16 N data 
are divided into zenith angle bins to show the relative 
asymmetry of the energy scale. These show asymmetries 
within ±0.5% (Figure]?]) and are similar with SK-I values. 

Quantitative representations of trigger efficiencies arc 
also obtained from 16 N data. The lowest threshold where 
the SK-II trigger is 100% efficient is 6.5 MeV whereas the 
SK-I threshold is 4.5 MeV. Figure [8] shows the trigger 
efficiencies for SK-II. A systematic error is assigned to 
the trigger efficiency by comparing the value given by 
16 N data and MC-simulated trigger events (0.5% on the 
total flux measurement). 



III. DATA ANALYSIS 

A. Trigger Scheme 

Like SK-I, SK-II has two levels of triggering for solar 
neutrino analysis: low energy [LE] and super-low energy 
[SLE] thresholds which require a minimum of about 14 
and 10 hit PMTs respectively to register an event. At the 
beginning of data taking in December 2002, only the LE 
trigger threshold was applied. At and above 8.0 MeV is 
where the LE trigger is 100% efficient. Later, the thresh- 
old was lowered and SLE data was taken with 100% ef- 
ficiency at 6.5 MeV. Ultimately, the LE+SLE analysis 
threshold was set to 7.0 MeV due to the large number of 



background events below this level. The LE analysis pe- 
riod lasted from December 24th, 2002 to July 14th, 2003 
for an exposure of 159 live days. The LE+SLE period 
lasted from July 15th, 2003 until October 5th, 2005 for 
an exposure of 632 live days. 

SLE triggered events are filtered online to reduce the 
amount of data written to limited storage space. Events 
reconstructed outside the fiducial volume are rejected. 
The data are reduced by a factor of approximately six. 
See Fig [2] for the vertex resolution as a function of energy 
of the online fitter. 



B. Background Reduction 

For SK-II, we implement a new two-part cut of events 
based on defined goodness functions of PMT timing and 
hit patterns. 

Many background events remain due to mis- 
reconstruction after the usual two-meter fiducial vol- 
ume cut, (which reduces background coming from the 
PMTs and blast shields.) Whereas a gamma-ray cut 
solely relies on vertex and directional reconstruction, the 
2-dimensional timing-pattern cut removes those events 
whose reconstruction should not be trusted. An opti- 
mized hit PMT timing goodness is defined (Equation 3.1l 
by comparing two timing residual Gaussian distributions, 
one with a width of a = 5 ns to encompass selected hits 
and the other with a to — 60 ns width characteristic of 
the PMT timing resolution for a single photo-electron: 



9t{v) 



Se 



(")- f O ^2 



(3.1) 



The effective hit time is defined as Ti{v) = U — \v — hi\/c 
which is just the timing residual t res id U ai of Equation |2.1| 
with added io- The sums are over all hits. 

The hit pattern goodness allows us to identify 
Cherenkov events by their azimuthal-symmetric ring pat- 
tern from a reconstructed vertex and direction. All others 
arc labeled mis-reconstructed or non-Cherenkov events. 
A goodness function g p (v) is defined for all directional 
events. 

A cut on the goodness values is made in tandem using 
a hyperbolic radius of g\ — > 0.25 and rejecting all 
other events. Figure [9] shows this background reduction 
cut on data and 8 B Monte Carlo in the 7.0-7.5 MeV bin. 
When the cut is applied to LINAC data and MC, a total 
flux systematic error of ±3.0% is obtained. 

The timing-hit pattern cut is treated as a second re- 
duction after the removal of noise and spallation events 
from the initial data set. All cuts and their efficiencies 



are shown in Figure 10 



6 



Si f f 




10 



ii 



12 14 
Energy (MeV) 



cn 

>■ 0.01 

I ° 

O -0.01 

(i -0.02 

— -0.03 
-0.04 



ui 0.015 
o 

5 o.oi 

o 

g. 0.005 
m 



8-0.005 
6 

s -0.01 

-0.015 



-II 

1(1) 


III 

(2; 


i i i i i i i i i i i i i i i i i i 

(3) 


1 


III. 
(5) 1 


- « . 

r ♦ 


* * 


+ 


♦ 

4 


+ 


+ 

it 


1 


+ 1 
4 : 


: i i 





♦ faj (6; fe) 




t 

+ 

1 


i i : 



(MOtMlDCMOO 



SOURCE z-position (m) 



l I I | 1 1 1 | I 1 1 | 1 1 I | I I 1 | I I I | 1 1 1 | I I I | 1 1 1 | I I 




t 


t 


T 


+ ; 


T 


, 1 , , , I , , , 




, i , , , i , 


, , i , , , : 



1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 

cos(zenith angle) 



FIG. 7: Top left: Deviation in energy scale between LINAC MC and data. Bottom left: differences in energy resolution between 
LINAC MC and data. Refer to Figurejfjjto relate data points with tank positions. Top right: 16 N energy scale deviation from 
MC. Representing the varying times calibration data were taken, (1) Nov. 2003, (2) March 2004, (3) July 2004, (4) Nov. 2004, 
(5) Sep. 2005. (a), (b), and (c) represent the x positions 15.20 m, 10.96 m, —14.49 m respectively. All other calibration data 
were taken at x = 0.35 m. Bottom right: 16 N energy scale deviation from MC for 6 selected zenith angles of the detector (-1 
is down). 




electron total energy (MeV) 

FIG. 8: Trigger efficiency as a function of energy. The black 
dots are 16 N calibration data and the lines are the best-fit 
error functions to the data (red solid is LE and green dashed 
is SLE). 




vertex (timing) goodness 



FIG. 9: PMT timing and hit pattern cut. Data (left) show 
an excess of mis-reconstructed and non-Chrerenkov events to 
the upper-left of the diagonal cut line. Approximately 78% 
(8%) of data (MC) events between 7.0-7.5 MeV are rejected 
by the cut. The color scale is to show the relative (normalized) 
number of events. 



C. Total Flux Result 



The SK-II solar neutrino signal is extracted from the 
strongly forward biased direction of recoil electrons from 
v-c elastic scattering. A likelihood fit to the signal and 



background is utilized to determine the flux. For a live 
time of 791 days of SK-II data from 7.0 to 20.0 MeV, the 
extracted number of signal events is 7212.8^5o'g(stat.) 
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FIG. 10: Summary of the data reduction steps (top) and their 
efficiencies on MC (bottom) 



TABLE I: SK-II systematic error of each item in %. Numbers 
in parentheses are the values obtained from calibration data 
before application to the neutrino flux. 
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FIG. 11: The angular distribution of solar neutrino candidate 
events. The flat line seen under the peak in the solar direction 
represents background contributions. 



D. Time Variation Results 

(2.38±0.05(stat.)iai5(sys.)) x 10 6 cm _2 sec _1 . 

1. Day-Night and Seasonal Variation 



It is statistically consistent with the SK-I value of (2.35± 
0.02(stat.) ±0.08(sys.)) x 10 6 cm^sec" 1 . The system- 
atic uncertainties of SK-I and II are mostly uncorrclatcd 
due to differences in energy scale, event selection, event 



reconstruction methods, etc. Figure 11 shows the angu- 
lar distribution of extracted solar neutrino events. Table 
I lists the SK-II systematic errors assigned for the total 
flux and day-night difference. 



Time variations of the solar neutrino flux are also de- 
termined by looking at day and night fluxes and the 
change in total flux at regular intervals during the live 
time of SK-II. The day and night fluxes are measured by 
selecting events which occur when the cosine of the solar 
zenith angle is less than zero (day) and greater than zero 
(night). Unlike the total flux, the day and night fluxes 
are quoted using a threshold of 7.5 MeV due to low sig- 
nal to noise ratio for the 7.0-7.5 MeV bin in the solar 
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direction after the data set is divided. Their values are 

<$> day = (2.31 ± 0.07(stat.) ± 0.15(sys.)) x 10 6 cm^sec" 
$night = (2.46±0.07(stat.)±0.16(sys.)) x 10 6 cm~ 2 sec" 

With these fluxes, the asymmetry value is found from 

A= day ~ ® night) /{\{® day + ® night))- The SK-II day- 
night difference yields 

A = -0.063 ± 0.042(stat.) ± 0.037(sys.). 

As with SK-I (A = -0.021 ±0.020(stat.)lg$ij(sys.)), no 
day-night asymmetry is discerned from the SK-II solar 
data set. The SK-I asymmetry value is statistically con- 
sistent with SK-II. 
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FIG. 12: Time dependence of the solar neutrino flux. The 
black points are from the 1496-day SK-I data set at a thresh- 
old of 5.0 MeV. The blue points are from the 791-day SK-II 
data set at a threshold of 7.0 MeV. The black line represents 
the expected 1/r 2 flux variations due to the eccentricity of 
the earth's orbit around the sun. Errors are statistical only. 
The absence of data points between SK-I and SK-II indicates 
dead time while construction of SK-II was occurring. 



the full period of solar cycle 23. To address any possible 
correlation of solar neutrino flux with sun spot number, 
,the SK-I and II flux time variation data are compiled in 
1-year bins between 1996 and 2006. The SK-I data set 
"(from 1996 to 2001) is taken from a 5.0 MeV threshold 
while SK-II is from 7.0 MeV. Errors are statistical only. 
From 1996 to the end of the SK-II phase in October 2005, 
the solar neutrino flux is stable and shows no pattern of 
correlation with the minima and maximum of solar cy- 
cle 23. This is consistent with (and a continuation of) 
the Kamiokande measurement and comparison with so- 
lar cycle 22 [TUj, albeit with a greater level of precision 
for Super-Kamiokande. 



» 3 

CM 

E 

o 

O 

£2.5 
x 

3 



1.5 K 



200 m 

H— » 

O 
Q. 

V) 

150 i 

(A 



1002 



50 



1998 2000 



2002 2004 2006 
YEAR 



FIG. 13: Time variation of the solar neutrino flux overlaid 
with sun spot number for solar cycle 23. Errors are statistical 
only. The SK-I and II 1-year binned solar flux data gives 
an agreement of \ 2 = 6.11 (52% c.l.) when compared to a 
straight line. 



E. Energy Spectrum 



The total flux variation as a function of time, or sea- 
sonal variation, for both SK-I and SK-II solar data is 



shown in Figure 12 Each bin represents 1.5 months and 
is seen to follow a sinusoidal trend consistent with the ex- 
pected 1 /r 2 flux variations due to the eccentricity of the 
earth's orbit around the sun. SK-II has excellent agree- 
ment with SK-I data, thus showing the continuation of 
the SK solar neutrino measurement through two phases 
of the detector. 



The recoil electron energy spectrum is obtained by di- 
viding the total flux into 17 energy bins ranging from 
7.0 to 20.0 MeV. The bin boundaries and flux values are 



listed in Table II. Figure 14 shows the observed energy 
spectrum divided by the expected spectrum without os- 
cillation determined from the BP2004 SSM [4\ The line 
through the spectrum represents the total SK-I 1496-day 
average. Again, as with the seasonal variation, SK-II 
shows excellent agreement with SK-I. 



2. Flux Correlation with Solar Activity 

With the completion of SK-II, the solar neutrino 
flux measurement of the Super-Kamiokande experiment 
spans an interval of 9.5 years. This closely coincides with 



IV. SK-II OSCILLATION ANALYSIS 

A. x 2 Minimization 

Oscillations of solar neutrinos have been studied by nu- 
merous experiments and have placed increasingly strin- 
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TABLE II: SK-II observed energy spectra expressed in units of event/kton/year. The errors in the observed rates are statistical 
only. The 7.0-7.5 MeV energy bin is excluded from the day-night analysis. Correction is made for the reduction efficiencies in 
Figure [TO] The expected rates neglecting oscillation are for the BP2004 SSM flux values. 8 Z is the angle between the z-axis of 
the detector and the vector from the Sun to the detector. 
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FIG. 14: Ratio of observed and expected energy spectra. The 
purple lines represent a ±1 sigma level of the energy correlated 
systematic errors. The black line represents the SK-I 1496- 
day average and shows agreement with SK-II. 



gent constraints on the mixing angle between neutrino 
mass and flavor eigenstates as well as neutrino mass dif- 
ference. In the statistically large data sample of SK-I, 
those constraints, assuming two flavor oscillations, favor 
the large mixing angle [LMA] region at 95% confidence 
level. The best fit values are given in the LMA region 
at tan 2 61=0.52 and Am 2 = 6.3 x lCT 5 eV 2 . The favored 
regions and corresponding best fit value are from a fit to 



the SK-I spectrum and time variation rates. The 8 B flux 
is also constrained by the SNO total rate [12] . 

The determination of the solar neutrino oscillation pa- 
rameters (#i2,Ami2) in SK-II is accomplished in much 
the same way as the previous SK-I result. Two neu- 
trino oscillation is assumed and for each set of oscillation 
parameters, a x 2 minimization of the total 8 B and hep 
neutrino flux is fit to the data. The entire SK-II ob- 
served spectrum is utilized from a 7.0 MeV threshold. 
The expected oscillated 8 B and hep flux is derived from 
numerically calculated MSW [5] v e survival probabilities 
and the unoscillated flux provided by the SSM. It is then 
converted to an expected SK-II rate spectrum by uti- 
lizing the v — e clastic scattering cross section and the 
SK-II detector's energy resolution. To account for the 
systematic uncertainties in energy resolution as well as 
the energy scale and the 8 B neutrino spectrum model 
shape, the combined rate predictions are modified by en- 
ergy shape factors, f{Ei,d B ,Ss,5 R ). 8 B , Ss, and 5 R rep- 
resent unfeertainty in the 8 B neutrino spectrum, SK-II 
energy scale, and SK-II energy resolution respectively. 
The function / serves to shift the rate predictions corre- 
sponding to a given uncertainty 5 in the data rate. The 
following equation shows the SK-II spectrum x 2 along 
with energy correlated systematic error shape factors ap- 
plied to the expected rate: 



2 _ {dj - (f3b z + rih t ) x f(Ej, S B , 5 S , S R )) 2 
Xsk-ii - / , a 2 

i=l 1 
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5r 
or 



■2Alog(£), (4.1) 



where di is the observed rate divided by the expected, 
unoscillated rate for the i th energy bin. Similarly, bi and 
hi are the predicted MSW oscillated rates divided by the 
unoscillated rate for 8 B and hep neutrinos respectively. 
(3 (77) scales the 8 B (hep) neutrino flux. The last term is 
the unbinned time-variation likelihood to the SK-II solar 
zenith angle flux variation above a 7.5 MeV threshold. 
This likelihood is analogous to the one used in SK-I. 

The energy uncorrelated systematic uncertainty is as- 
signed to 4.8% (the quadrature sum of the energy inde- 
pendent errors in Table I) and is used to describe the 
error on the total rate. For the spectrum rate uncertain- 
ties, the value 4.8% is conservatively assigned to each bin 
and is added in quadrature to the statistical error (Table 
II) to equal Oi in the SK-II % 2 . See the appendix of jTTj 
for more details. 



B. Oscillation Results - SK Only 

A minimization of the x 2 m the previous section yields 
excluded regions when (3 and 77 are left unconstrained. 
By constraining the 8 B flux to the total NC flux value 
from SNO [T2], allowed parameter regions can be ob- 



tained. Figure 15 shows both excluded and allowed re- 
gions at 95% confidence level. They are consistent with 
previous SK-I results. The primary constraint in SK-II 
is from the time-variation data although some spectral 
exclusion is also seen at Am 2 « 10 _4 eV 2 . The same os- 
cillation analysis is performed while including \ 2 terms 
corresponding to the SK-I values (namely, the spectrum 
and unbinned time variation for SK-I). In this combined 
analysis, SK-II helps expand the 95% c.l. exclusion from 
the SK-I-only analysis, mostly along a region dominated 
by the spectral constraint (10 -4 < tan 2 6 < 0.4 and 
4 x 10~ 5 eV 2 < Am 2 < 2 x 10" 4 eV 2 ). However, when 
constraining 8 B to the SNO NC flux, the SK allowed 
regions are largely unaffected by the addition of SK-II 
data. 



C. Oscillation Results - SK and Other Solar 
Experiments 

The combination of other solar neutrino experiments 
such as the SNO and radiochemical results with the 



SK combined analysis is accomplished by fitting the to- 
tal CC and NC rates observed by SNO's 306-day pure 
D 2 [T3] and 391-day salt phases [IS]. Also, the SNO 
NC constrained predicted day-night asymmetry for the 
pure D 2 phase is used for added exclusion power. The 
radiochemical experiments of Homestake, GALLEX, and 
SAGE [T3] are then added using the best 8 B and hep 
fluxes from the SK-SNO fit. Figure [T5] shows the com- 
bined solar allowed areas. The best fit parameter set is 
tan 2 = 0.40 and Am 2 = 6.03 x 10~ 5 eV 2 consistent with 
the SK-I global analysis. 



V. CONCLUSION 

Super-Kamiokande has measured the solar 8 B flux to 
be (2.38±0.05(stat.)tjj;^(sys.)) x 10 6 cm^sec" 1 during 
its second phase. The uncertainties in SK-II are larger 
than in SK-I but a low analysis threshold of 7 MeV was 
achieved (7.5 MeV in the day-night variation analysis). A 
day-night asymmetry value was observed to be — 0.063 ± 
0.042(stat.) ± 0.037(sys.) which is consistent with zero 
and the result from SK-I. SK-II has brought the total SK 
time-dependent flux measurement to a length of 9.5 years 
and this measurement is compared with solar activity in 
solar cycle 23 resulting in no strong correlation. In the 
combined SK-I and SK-II global oscillation analysis, the 
best fit is found to favor the LMA region at tan 2 9 = 
0.40 and Am 2 = 6.03 x 10~ 5 eV 2 , in excellent agreement 
with previous solar neutrino oscillation measurements. 
SK-I and SK-II agree well, showing no evidence of any 
systematic effects from the introduction of new methods, 
blast shields, reduced PMT coverage, etc. 
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